Phenol Molecular Sheets Woven by Water Cavities in Hydrophobic Slit Nanospaces.

Despite extensive research over the last several decades, the microscopic characterization of topological phases of adsorbed phenol from aqueous solutions in carbon micropores (pore size < 2.0 nm), which are believed to exhibit a solid and quasi-solid character, has not been reported. Here, we present a combined experimental and molecular level study of phenol adsorption from neutral water solutions in graphitic carbon micropores. Theoretical and experimental results show high adsorption of phenol and negligible coadsorption of water in hydrophobic graphitic micropores (super-sieving effect). Graphic processing unit-accelerated molecular dynamics simulation of phenol adsorption from water solutions in a realistic model of carbon micropores reveal the formation of two-dimensional phenol crystals with a peculiar pattern of hydrophilic-hydrophobic stripes in 0.8 nm supermicropores. In wider micropores, disordered phenol assemblies with water clusters, linear chains, and cavities of various sizes are found. The highest surface density of phenol is computed in 1.8 nm supermicropores. The percolating water cluster spanning the entire pore space is found in 2.0 nm supermicropores. Our findings open the door for the design of better materials for purification of aqueous solutions from nonelectrolyte micropollution.

Nanoscale Water Contact Angle on Polytetrafluoroethylene Surfaces Characterized by Molecular Dynamics-Atomic Force Microscopy Imaging.

The aim of this study is to link polytetrafluoroethylene (PTFE) surface characteristics with its wetting properties in the nanoscale. To do this using molecular dynamics (MD) simulation, three series of rough PTFE surfaces were generated by annealing and compressing and next characterized by the application of the MD version of the atomic force microscopy (AFM) method. The values of specific surface areas were additionally calculated. The TIP4P/2005 water model was used to study the wetting properties of obtained PTFE samples. The simulated water contact angle (WCA) value for the most flat (but slightly rough) sample having PTFE density is equal to 106.94°, and it is close to the value suggested for a perfect PTFE surface on the basis of experimental results. Also, the changes in the WCA with PTFE compression are in the same range as experimentally reported. The obtained MD simulation results make it possible to link, for the first time, the WCA values with the surface MD-AFM root-mean-square roughness and with the PTFE density. Finally, we show that for PTFE wetting in the nanoscale, the line tension is negligible and the Bormashenko's equation reduces to the Cassie-Baxter (CB) model. In fact, our simulation results are close to the CB mechanism.

Nanoscale Insight into the Mechanism of a Highly Oriented Pyrolytic Graphite Edge Surface Wetting by "Interferencing" Water.

The new molecular dynamics simulation results showing the influence of the edge carbon surface atoms on the wettability of a highly oriented pyrolytic graphite (HOPG) surface with water nanodroplets are reported. The conditions for the occurrence of the Wenzel effect are discussed, and the Cassie-to-Wenzel transition (CTWT) mechanism in the nanoscale is explored. This transition is detected by the application of a new procedure showing that the CTWT point shifts toward larger values of carbon-oxygen potential well depth with the decrease in the HOPG side angle. It is concluded that the Wenzel effect significantly contributes to the contact angles (CAs) measured for the HOPG surfaces. The Wenzel effect is also very important for the "HOPG" structures possessing the disturbed C-C interlayer distance, and its influence on the water nanodroplet CAs is strongly pronounced. The structure of water confined inside slits and on a HOPG surface is studied using the analysis of the density profiles, the number of hydrogen bonds, and, modified for the purpose of this study, structure factor. The detailed analysis of all parameters describing confined water leads to the conclusion about the presence of characteristic interference patterns revealed as a result of long-term simulation. A simple model describing this effect is proposed as the starting point for further considerations.

Water Adsorption Property of Hierarchically Nanoporous Detonation Nanodiamonds.

The detonation nanodiamonds form the aggregate having interparticle voids, giving a marked hygroscopic property. As the relationship between pore structure and water adsorption of aggregated nanodiamonds is not well understood yet, adsorption isotherms of N2 at 77 K and of water vapor at 298 K of the well-characterized aggregated nanodiamonds were measured. HR-TEM and X-ray diffraction showed that the nanodiamonds were highly crystalline and their average crystallite size was 4.5 nm. The presence of the graphitic layers on the nanodiamond particle surface was confirmed by the EELS examination. The pore size distribution analysis showed that nanodiamonds had a few ultramicropores with predominant mesopores of 4.5 nm in average size. The water vapor adsorption isotherm of IUPAC Type V indicates the hydrophobicity of the nanodiamond aggregates, with the presence of hydrophilic sites. Then the hygroscopic nature of nanodiamonds should be associated with the surface functionalities of the graphitic shell and the ultramicropores on the mesopore walls.

Water nanodroplet on a graphene surface-a new old system.

The major subject of our study is the accuracy of contact angle calculations. Reporting new simulation data for graphene-water systems, we show that the majority of previously reported data should be treated with caution, since the proper contact angle can be recorded only after a sufficiently long simulation time. It has been proven that-if one wants to gain accuracy greater than 0.1°-long calculations (exceeding 50 ns) are required. Finally, we also show, using both a Groningen Machine for Chemical Simulations (GROMACS) package and our new molecular dynamics (MD) code, that the changes in the contact angle, caused by graphene bottom layer rotation, are within the range of calculation error. We also propose a novel definition of the bottom of the droplet as the height where the density is half the density of liquid water. This new definition is applied in the method of the contact angle calculation from the MD simulation data.

To what extent can mutual shifting of folded carbonaceous walls in slit-like pores affect their adsorption properties?

We have performed systematic Monte Carlo studies on the influence of shifting the walls in slit-like systems constructed from folded graphene sheets on their adsorption properties. Specifically, we have analysed the effect on the mechanism of argon adsorption (T = 87 K) and on adsorption and separation of three binary gas mixtures: CO2/N2, CO2/CH4 and CH4/N2 (T = 298 K). The effects of the changes in interlayer distance were also determined. We show that folding of the walls significantly improves the adsorption and separation properties in comparison to ideal slit-like systems. Moreover, we demonstrate that mutual shift of sheets (for small interlayer distances) causes the appearance of small pores between opposite bulges. This causes an increase in vapour adsorption at low pressures. Due to overlapping of interactions with opposite walls causing an increase in adsorption energy, the mutual shift of sheets is also connected with the rise in efficiency of mixture separation. The effects connected with sheet orientation vanish as the interlayer distance increases.

MD simulation of organics adsorption from aqueous solution in carbon slit-like pores. Foundations of the pore blocking effect.

The results of systematic studies of organics adsorption from aqueous solutions (at the neutral pH level) in a system of slit-like carbon pores having different sizes and oxygen groups located at the pore mouth are reported. Using molecular dynamics simulations (GROMACS package) the properties of adsorbent-adsorbate (benzene, phenol or paracetamol) as well as adsorbent-water systems are discussed. After the introduction of surface oxygen functionalities, adsorption of organic compounds decreases (in accordance with experimental data) and this is caused by the accumulation of water molecules at pore entrances. The pore blocking effect decreases with the diameter of slits and practically vanishes for widths larger than approx. 0.68 nm. We observed the increase in phenol adsorption with the rise in temperature. Moreover, adsorbed molecules occupy the external surface of the slit pores (the entrances) in the case of oxidized adsorbents. Among the studied molecules benzene, phenol and paracetamol prefer an almost flat orientation and with the rise in the pore width the number of molecules oriented in parallel decreases. The decrease or increase in temperature (with respect to 298 K) leads to insignificant changes of angular orientation of adsorbed molecules.

Testing isotherm models and recovering empirical relationships for adsorption in microporous carbons using virtual carbon models and grand canonical Monte Carlo simulations.

Using the plausible model of activated carbon proposed by Harris and co-workers and grand canonical Monte Carlo simulations, we study the applicability of standard methods for describing adsorption data on microporous carbons widely used in adsorption science. Two carbon structures are studied, one with a small distribution of micropores in the range up to 1 nm, and the other with micropores covering a wide range of porosity. For both structures, adsorption isotherms of noble gases (from Ne to Xe), carbon tetrachloride and benzene are simulated. The data obtained are considered in terms of Dubinin-Radushkevich plots. Moreover, for benzene and carbon tetrachloride the temperature invariance of the characteristic curve is also studied. We show that using simulated data some empirical relationships obtained from experiment can be successfully recovered. Next we test the applicability of Dubinin's related models including the Dubinin-Izotova, Dubinin-Radushkevich-Stoeckli, and Jaroniec-Choma equations. The results obtained demonstrate the limits and applications of the models studied in the field of carbon porosity characterization.

How realistic is the pore size distribution calculated from adsorption isotherms if activated carbon is composed of fullerene-like fragments?

A plausible model for the structure of non-graphitizing carbon is one which consists of curved, fullerene-like fragments grouped together in a random arrangement. Although this model was proposed several years ago, there have been no attempts to calculate the properties of such a structure. Here, we determine the density, pore size distribution and adsorption properties of a model porous carbon constructed from fullerene-like elements. Using the method proposed recently by Bhattacharya and Gubbins (BG), which was tested in this study for ideal and defective carbon slits, the pore size distributions (PSDs) of the initial model and two related carbon models are calculated. The obtained PSD curves show that two structures are micro-mesoporous (with different ratio of micro/mesopores) and the third is strictly microporous. Using the grand canonical Monte Carlo (GCMC) method, adsorption isotherms of Ar (87 K) are simulated for all the structures. Finally PSD curves are calculated using the Horvath-Kawazoe, non-local density functional theory (NLDFT), Nguyen and Do, and Barrett-Joyner-Halenda (BJH) approaches, and compared with those predicted by the BG method. This is the first study in which different methods of calculation of PSDs for carbons from adsorption data can be really verified, since absolute (i.e. true) PSDs are obtained using the BG method. This is also the first study reporting the results of computer simulations of adsorption on fullerene-like carbon models.

Hyper-parallel tempering Monte Carlo simulations of Ar adsorption in new models of microporous non-graphitizing activated carbon: effect of microporosity.

The adsorption of gases on microporous carbons is still poorly understood, partly because the structure of these carbons is not well known. Here, a model of microporous carbons based on fullerene-like fragments is used as the basis for a theoretical study of Ar adsorption on carbon. First, a simulation box was constructed, containing a plausible arrangement of carbon fragments. Next, using a new Monte Carlo simulation algorithm, two types of carbon fragments were gradually placed into the initial structure to increase its microporosity. Thirty six different microporous carbon structures were generated in this way. Using the method proposed recently by Bhattacharya and Gubbins (BG), the micropore size distributions of the obtained carbon models and the average micropore diameters were calculated. For ten chosen structures, Ar adsorption isotherms (87 K) were simulated via the hyper-parallel tempering Monte Carlo simulation method. The isotherms obtained in this way were described by widely applied methods of microporous carbon characterisation, i.e. Nguyen and Do, Horvath-Kawazoe, high-resolution α(s) plots, adsorption potential distributions and the Dubinin-Astakhov (DA) equation. From simulated isotherms described by the DA equation, the average micropore diameters were calculated using empirical relationships proposed by different authors and they were compared with those from the BG method.